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Results of DC and frequency dependent conductivity in the quantum limit, i.e. Tiu) > ksT, for a 
broad range of dopant concentrations in nominally uncompensated, crystalline phosphorous doped 
silicon and amorphous niobium-silicon alloys are reported. These materials fall under the general 
category of disordered insulating systems, which are referred to as electron glasses. Using microwave 
resonant cavities and quasi-optical millimeter wave spectroscopy we are able to study the frequency 
dependent response on the insulating side of the metal-insulator transition. We identify a quantum 
critical regime, a Fermi glass regime and a Coulomb glass regime. Our phenomenological results 
lead to a phase diagram description, or taxonomy, of the electrodynamic response of electron glass 
systems. 



PACS numbers: 72.20.Ee, 71.30.+h, 71. 45. Cm 



I. INTRODUCTION 

Disordered insulating systems have garnered much 
interest from theorists and experimentalists alike over 
many years. A plethora of results exist on the low energy 
physics of these disordered insulating systems from DC 
conductivity experiments, but it is only recently that in- 
vestigations of the T — > frequency dependent response 
have been performed even though theoretical predictions 
have existed for many years. We have measured the 
DC conductivity as well as the complex frequency de- 
pendent conductivity, the latter being measured in the 
quantum- limit, i.e. Tito > ksT for a few different disor- 
dered insulating systems. Measurements were performed 
on materials that fall on the insulating side of the metal 
insulator transition (MIT), with a major focus on the 
effects of electron-electron interactions. Theoretical pre- 
dictions concerning the charge transport in systems with 
carriers localized in the Anderson sense exist both for sys- 
tems without electron-electron interactions, these mate- 
rials being referred to as Fermi glasses 1 , as well as for sys- 
tems that include the effects of the Coulomb interaction, 
these materials being referred to as Coulomb glasses 2 . 
The majority dopant concentration of doped crystalline 
semiconductors and amorphous metal-semiconductor al- 
loys can easily be varied using standard semiconductor 
growth techniques, and hence these materials provide a 
perfect playground to study the low energy electrody- 
namics of localized charge carriers. 

We focus our investigation on a broad region of phase 
space from close the the quantum critical point, i.e. the 
critical dopant concentration at which the MIT occurs, 
x c , to deep within the insulating regime. Recent work 
on the AC conductivity of Si:B 3 and Si:P 4 have shown a 
crossover from the Fermi glass to the Coulomb glass-like 
behavior qualitatively consistent with the theoretical ex- 
pectations, but in all cases, the crossover is much sharper 
than predicted by theory. There exist two relevant en- 
ergy scales in this problem. The first is the Coulomb 



interaction energy, U, which Efros and Shklovskii (ES) 
predicted determines the crossover energy between char- 
acteristic Fermi glass and Coulomb glass- like behavior 2 . 
An important result of the interactions is a depletion 
in the density of states of the charge carriers about the 
Fermi level termed a Coulomb gap A. In a recent result, 
Lee et al. claim from their measurements on Si:B that 
the crossover energy scale is determined by the Coulomb 
gap width, A, as opposed to the theoretically predicted 
Coulomb interaction energy, U. Arguments based on the 
concentration dependence of these two parameters ob- 
tained from measurements over a broad range of dopant 
concentrations in another canonical electron glass sys- 
tem, Si:P, have recently been published by us casting 
doubt on this interpretation 4 . This highlights the use- 
fulness and importance of investigating a wide range of 
majority dopant concentrations to fully understand these 
disordered insulating systems. 

The primary goals of this work then are to present 
new data and to draw on results of previous work on 
other electron glass systems 4-6 to look at the insulating 
side of the MIT using AC conductivity data and a broad 
range of dopant concentrations in order to come up with 
a phcnomenologically supported phase diagram outlining 
the general classification scheme for the electrodynamic 
response of electron glasses. DC conductivity data are 
presented as well and a side by side comparison of the AC 
and DC data provides for unique insight. Furthermore 
the AC data allows for a detailed comparison of the elec- 
tron glass regime with theoretical predictions. We find 
that some of the predicted features are observed, such 
as the frequency dependent conductivity power laws. In 
fact it is these powers that provide the taxonomic infor- 
mation in the phase diagram to be presented herein. We 
also note discrepancies that arise in all the reported mea- 
surements, such as the sharper than predicted crossovers 
from Fermi glass to Coulomb glass behavior. 
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II. BACKGROUND 
A. The MIT as a Quantum Phase Transition 

The metal insulator transition is an archetypical quan- 
tum phase transition (QPT). The insulating portion of 
phase space associated with Fermi glass and Coulomb 
glass behavior borders the quantum critical regime. The 
quantum critical region of phase space is that region 
where quantum mechanical effects due to being in prox- 
imity to the zero temperature quantum critical point 
manifest themselves even at finite temperatures or fre- 
quencies. Specifically the electrodynamic response in a 
disordered electronic system, such as the lightly doped 
insulating semiconductors studied in this work, should 
smoothly transform from Fermi glass or Coulomb glass- 
like behavior into the appropriate quantum critical dy- 
namic response as the dopant concentration is tuned and 
the MIT is approached. 

For the MIT, the dopant concentration which is di- 
rectly related to the charge carrier concentration, x, is 
typically the tuning parameter that determines whether 
at zero temperature the material is a metal, i.e. has 
a finite DC conductivity, or an insulator, i.e. the DC 
conductivity is zero. There exists necessarily then a 
critical concentration, x c , at which the phase transition 
occurs thereby defining a boundary between insulating 
and metallic ground states. The disorder induced metal- 
insulator transition and the dynamical behavior near the 
critical point has been extensively theoretically studied 
over the past years both with and without electron- 
electron interactions 7-10 . A thorough review has been 
written by Belitz and Kirkpatrick 11 . 

The condition for observing quantum critical behavior 
versus electron glass-like behavior on the insulating side 
of the MIT is determined by the relative size of certain 
characteristic length scales. The characteristic thermal 
length scale is 

It oc (ksT)- 1 /* (1) 

and the frequency dependent length scale is 

t u oc (Huj)- 1 /* (2) 

where z is the dynamical scaling exponent. These lengths 
are termed "dephasing" lengths by Sondhi et al. 12 . The 
localization length, a measure of how far the electronic 
wavefunction extends, is the other important length 
scale. It diverges as the transition is approached and 
can be written in terms of concentration as 

£oc \x-x c \- v (3) 

where v is the localization length exponent. 

At finite temperature or frequency an intermediate 
regime exists between the metallic and insulating regimes 
that displays quantum critical (QC) behavior as shown 



in Figure 1. In the QC regime the localization length 
is greater than the dephasing length scales, £t,u> < £• 
At the crossover from QC to the insulating regime this 
can be thought of heuristically as the electron scattering 
within a length scale that is shorter than the width of 
the wavefunction envelope that defines how localized that 
electron is. In other words the electron can undergo scat- 
tering akin to normal transport before " realizing" that it 
is localized. On the metallic side, the relevant length 
scale is commonly taken to be the screening length. In 
the opposite limit, £ < £t,u the system exhibits metal- 
lic or insulating behavior. For instance in the insulating 
regime, the material will display the appropriate hop- 
ping conduction consistent with the DC conductivity or 
the photon assisted frequency dependent conductivity. A 
schematic representation of the phase diagram showing 
the regime boundaries is shown in Figure 1. 
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FIG. 1. A schematic of the crossover diagram close to a 
quantum critical point showing the regime boundaries be- 
tween quantum critical, metallic and insulating behavior. The 
horizontal axis represents increasing concentration and the 
vertical axis represents either the temperature or frequency 
energy scale. The crossover from QC to either metallic or 
insulating behavior is determined by when the length scales 
are equivalent, i.e. £t,^ ~ f . 

From equations (1), (2) and (3) we can determine 
a functional form, in terms of concentration, for the 
crossover between the QC and the metallic or insulat- 
ing regimes by setting the dephasing length equal to the 
relevant length scale, £t,u ~ £• In so doing we find 

hujQc <x\x- x c \ zu . (4) 

In this work the approach to criticality from the insulat- 
ing side of the MIT for two types of disordered insulators, 
Si:P and NbSi, will be compared and discussed. 
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B. DC Variable Range Hopping and 
the Coulomb Gap 

Mott first noted that at low temperatures in a lightly 
doped disordered semiconductor a form of DC con- 
duction unique from nearest-neighbor hopping would 
occur 13 . It is commonly referred to as variable range hop- 
ping (VRH). Mott's derivation did not include electron- 
electron interactions and thus Mott-typc hopping de- 
scribes the variable range hopping mechanism in a non- 
interacting Fermi glass 1 . In Mott's theory: 

a^ozexpi-iTo/T) 1 /^] (5) 

where d is the effective spatial dimension for the hopping 
electrons. The characteristic Mott temperature is given 

by 

To = f3/k B N a Z d . (6) 

It is the characteristic energy scale of the level spacing 
in a volume £ d , and dependent on the density of states 
(DOS), N . Here (5 is a constant of order one. 

The above derivation was based on the assumption 
of an approximately constant density of states near 
the Fermi level, but as pointed out by Pollak 14 and 
Ambegaokar 15 , this assumption is invalid when corre- 
lation effects are included. Upon including electron- 
electron interactions one finds that the Coulombic repul- 
sion is responsible for opening up a gap at the Fermi en- 
ergy, and that the functional form for the DC conductiv- 
ity takes on an alternate temperature dependence. This 
gap, termed the Coulomb gap by Efros and Shklovskii 16 
is given by 



where S\ is the full dielectric constant. The Coulomb gap 
is a "soft" gap in the density of states near the Fermi 
level created by the Coulomb interaction between local- 
ized electrons. (The term soft is used because the density 
of states vanishes only at the Fermi energy.) 

Efros and Shklovskii followed Mott's derivation for the 
DC temperature dependent VRH conductivity, but in- 
stead of using a constant density of states around the 
Fermi level, they incorporated the Coulomb gap and 
therefore necessarily took correlation effects into account. 
The appellation Coulomb glass describes an electron 
glass system which incorporates Coulombic interactions. 
In so doing ES found the following form for the DC tem- 
perature dependent VRH conductivity 

a dc cx expHU/T) 1 / 2 } (8) 

where the characteristic ES temperature is given by 

= 2.8e 2 /k B e^. (9) 



The prefactor 2.8 is calculated from percolation threshold 
models 17 ' 18 . Note that the power in Eq. (8) is i in all 
dimensions unlike in Mott's formalism. 

For the systems we are studying, any given sample with 
a particular dopant concentration may possibly exhibit 
more than one type of behavior. Depending on the ener- 
gies at which one is probing the system one might observe 
behavior consistent with Fermi glass, Coulomb glass or 
QC-like behavior. Therefore, by investigating a broad 
range of frequency or temperature space, our experimen- 
tal window might be broad enough to see crossovers from 
one type of behavior to another. 

The most indicative piece of experimental evidence 
that most authors cite as to whether the materials be- 
ing studied are Fermi glasses exhibiting Mott VRH or 
Coulomb glasses exhibiting ES VRH is the DC conduc- 
tivity. The vicissitudes of using this experimental data 
to claim a material is a Coulomb or Fermi glass lie in 
the fact that only below a certain temperature does a 
material enter the VRH regime, whereas above this cut- 
off temperature other forms of activated hopping occur. 
Thus only a small portion of low temperature data will 
show any semblance of VRH behavior, and in three di- 
mensions, a small section of temperature dependent data 
may fit a power law of \ or j almost equally as well mak- 
ing it difficult to say with any certainty which mechanism 
of charge excitation is occurring. Measurements of the 
T — > K frequency dependent conductivity then allows 
for new insight into the MIT and the electrodynamics of 
electron glasses. Measurements of the DC conductivity 
on the same samples as those on which the AC conductiv- 
ity measurements were taken offers a consistency check 
with previous investigations and allows for comparison of 
salient features. 



C. AC Photon-assisted Hopping 

Mott proposed a form for the T = K photon as- 
sisted frequency dependent hopping conductivity for the 
case of non-interacting electrons based on the theory of 
resonant absorption using a simple one electron model 
of a disordered system, commonly referred to as the pair 
approximation 1,19 . The resulting real part of the AC 
conductivity, <7i, in Mott's formalism is as follows: 

ai = ae 2 Nfe 5 hu 2 [ln(2I Q /huj)} 4 . (10) 

Here a is some constant close to one, and Iq is the pre- 
factor of the overlap integral (commonly taken to be the 
Bohr energy of the dopant 20 ). Most notably we see that 
in Mott's formalism, the conductivity has an approxi- 
mately quadratic frequency dependence, plus logarithmic 
corrections. 

By taking into account the mean Coulomb interaction 
between two sites forming a resonant pair 

U(r u ) = e 2 / El r u (11) 
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where r w = £[ln(2Io/fuj)} is the most probable hop dis- 
tance between pairs, Efros and Shklovskii derived the real 
part of the AC conductivity to be: 

ai = ae 2 A^ 5 ^/n(2VM] 4 [ftu> + U(r w )}. (12) 

This formula takes on a different frequency dependence in 
two limiting cases. In the case where the photon energy 
is greater than the Coulomb interaction energy, huj > 
U(r u ), one retrieves the same approximately quadratic 
frequency dependence in Mott's formula (10), and here 
the Coulomb glass is indistinguishable from the Fermi 
glass in so far as the frequency dependent conductivity 
is concerned. For the alternate case i.e. when the pho- 
ton energy is smaller than the Coulomb interaction en- 
ergy, hcu < U(r u] ) 1 the conductivity of a Coulomb glass 
will show an approximately linear dependence on fre- 
quency, plus logarithmic corrections. Model data for this 
functional form is shown in Figure 2 showing that the 
crossover is predicted to occur over many decades of fre- 
quency. 




Frequency 



FIG. 2. A schematic of the crossover function from Mott 
to ES-type AC conduction set forth by Efros and Shklovskii 
and defined in Eq. (12) displayed on logarithmic axes. The 
vertical axis is conductivity with arbitrary units and the hori- 
zontal axis has been set to show the crossover occurring in the 
millimeter-wave range if the units are taken to be GHz. The 
dashed and dotted straight lines show the asymptotically ap- 
proached linear and quadratic dependencies of the crossover 
function. 

As for the DC case, one further ramification of includ- 
ing electron-electron interactions is the depletion of the 
single particle density of states around the Fermi level, 
i.e. the Coulomb Gap A. In the special case when 
A > U(r u ) > Tied Efros and Shklovskii calculated the 
real part of the frequency dependent conductivity and 
obtained 

1 LUEl 

ai ~ lb~Zn(2/ /M' 

so that the dependence is slightly super-linear, assuming 
a frequency-independent dielectric constant, and inde- 



pendent of the localization length, which would indicate 
a very weak concentration dependence. 

Thus, by using the frequency dependent conductivity 
one has a unique means to determine the relevant energy 
scales between interacting and non-interacting behav- 
ior. Several recent experiments 3 ' 21 have shown a possible 
crossover from interacting to non-interacting behavior in 
the same material. Some of this experimental data has 
been interpreted to indicate that the Coulomb gap width 
energy scale, certainly another germane energy scale in 
the system, sets the crossover energy scale in opposition 
to the theoretical prediction that the Coulomb interac- 
tion energy dictates the crossover from Fermi glass to 
Coulomb glass-like behavior. These arguments though 
were based on measurements covering a limited range of 
dopant concentrations. Therefore the frequency depen- 
dent conductivity measured in the low energy excitation 
range for a broad range of concentrations can provide 
useful information about the crossover energy scale, be 
it the Coulomb gap width or the Coulomb interaction 
energy. 

III. EXPERIMENTS 

Experiments were performed on nominally uncompen- 
sated phosphorous doped silicon samples that were cut 
from a boule grown by Recticon Enterprises Inc. to a 
specification of 5 cm in diameter with a majority dopant 
gradient along the axis. This boule, grown using the 
Czochralski method, was subsequently sliced into 1 mm 
thick wafers. By the nature of the Czochralski method for 
growing a crystal, there is an inherent gradient in the ma- 
jority dopant concentration along the axis of the boule. 
A greater range of variation in this dopant concentration 
along the axis of the boule can be achieved by continu- 
ously adding progressively more dopant of choice, in our 
case phosphorous, to the molten silicon as the ingot is 
being drawn from the liquid mixture. As summarized in 
the first column in Table I, we obtained samples spanning 
a majority dopant concentration from 1.37 x 10 18 cm -3 
to 2.98 x 10 18 cmT z . The middle column lists the ratios 
of the samples as a percentage of the dopant concentra- 
tion in each sample relative to the critical concentration 
at which the MIT occurs, x c = 3.5 x 10 18 cm~ 3 . The fi- 
nal column lists the room temperature resistivity of each 
sample in (Q cm) measured using an ADE 6035 resistiv- 
ity gauge, which can be directly related to the dopant 
concentration, knowing the ilk of majority dopant in the 
silicon via a phenomenological scale set forth by Thurber 
ct al. 22 . 

Previously reported results from another electron glass 
system, namely amorphous NbSi, will be drawn upon 
when discussing the general features of the electron glass 
phase diagram. These samples were created by sputter- 
ing onto sapphire substrates, and both millimeter-wave 
and DC conductivity measurements were performed, 
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x xlO 18 cm,- 3 


Percent of Critical: — 


p dc (300 K) (Q cm) 


2.98 


85% 


0.0133 


2.66 


76% 


0.0142 


2.42 


69% 


0.0150 


2.17 


62% 


0.0159 


1.96 


56% 


0.0168 


1.74 


50% 


0.0179 


1.56 


45% 


0.0189 


1.37 


39% 


0.0202 



TABLE I. The material parameters for eight Si:P samples 
obtained from a boule grown by Recticon Enterprises Inc. 
The left hand column lists the concentration of the major- 
ity dopant phosphorous in each sample. The middle column 
lists the ratio as a percentage of the dopant concentration in 
each sample relative to the critical concentration of the MIT, 
namely x c — 3.5 x 10 18 cm -3 . The right hand column lists 
the value of the room temperature DC resistivity value for 
each sample. 

analogous to those done on Si:P. More detailed sample 
specifications for these amorphous metal-semiconductor 
alloys and experimental results can be found in the 
literature 5 . 

It is important to perform DC conductivity measure- 
ments on the same Si:P samples for which the AC mea- 
surements were done for a number of reasons. First, the 
DC measurements offer a consistency check to previous 
experimental results. Also, one can compare important 
parameters, such as the localization length, as deter- 
mined separately from the DC and AC data. These mea- 
surements also provided a low temperature value of the 
DC conductivity to which one can normalize the real part 
of the microwave resonant cavity complex conductivity. 
For example, at approximately 25 K the real part of the 
conductivity measured at the center frequencies of the 
resonant cavities, namely 35 GHz and 60 GHz, should 
be dominated by the thermally driven charge transport 
mechanism that determines the DC conductivity. 

DC conductivity measurements were performed on the 
Si:P samples using a standard four probe measurement 
technique from liquid 4 He temperature up to approxi- 
mately 100 K. Figure 3 shows data from three samples 
from 4.2 K up to 20 K. At the higher temperatures shown 
in the main portion of the figure for the the 69 % sam- 
ple, and also evident in the inset plot of the 56 % sample 
measured up to approximately 200 K, the curves show a 
rounding off indicative of a conduction process consistent 
with a Boltzmann type activation. Also evident in the 
middle of the displayed temperature range is a kink in 
each of the sample's data indicative of a crossover from 
one form of activation to another. Visually the crossover 
temperature, T x , for each of the samples decreases with 
increasing dopant concentration. This would be consis- 
tent with a crossover from Mott to ES-type VRH, but 
other experimental evidence, to be discussed presently, 
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FIG. 3. The DC conductivity versus temperature on loga- 
rithmic scales for a series of Si:P samples. At high tempera- 
tures one can notice a crossover to activated Boltzmann-type 
behavior. Also noticeable is a crossover temperature, T x , in 
the middle of the displayed temperature range that decreases 
with increasing dopant concentration. 

indicates that Mott and ES-type VRH theories do not 
correctly describe the upper and lower pieces of the DC 
conductivity data shown in Figure 3. 

Figure 4 shows the measured DC conductivity from 
the 62% sample both prior to and after etching in a 
4% HF + 96% HN0 3 solution for 20 minutes. Unless 
otherwise specified all AC and DC conductivity mea- 
surements reported herein were performed on etched Si:P 
samples. As first pointed out by Ootuka et al. 23 , a crys- 
talline semiconductor sample sliced from a boule using a 
diamond saw will create surface defects. As indicated 
in the figure, at low temperatures the surface defects 
seem to provide a channel for conduction with less resis- 
tance than the surface defect-free bulk sample, increasing 
the overall measured bulk conductivity. DC conductiv- 
ity measurements on samples etched for 20 minutes and 2 
hours showed no difference indicating that 20 minutes of 
etching sufficed to remedy the effects of surface defects, 
consistent with the findings of Ootuka et al. 23 . Regard- 
less of the presence of surface defects altering the relative 
magnitude of the measured bulk conductivity, the kink 
in the data indicating a crossover temperature, T x , is still 
clearly present in both sets of data. 

Both components of the complex frequency dependent 
conductivity were determined from transmission mea- 
surements using a technique called quasi-optical spec- 
troscopy. Quasi-optical spectroscopy has been utilized 
and partially developed by our research group to inves- 
tigate the electrodynamic response of materials in the 
millimeter to sub-millimeter wave range (i.e. the terra- 
hertz frequency regime) 24 ' 25 . Historically the millimeter 
and sub-millimeter wave spectral range, higher in fre- 
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FIG. 4. The DC conductivity versus temperature on a log- 
arithmic scale for the 62% Si:P sample both prior to and after 
etching the sample. 



quency than standard radio-frequency (RF) techniques 
and lower in frequency than common infrared (IR) tech- 
niques, has not been studied, leaving a gap in the knowl- 
edge of the optical response of physically intriguing sys- 
tems. This spectral range corresponds with the range of 
energies wherein electron-electron interactions are pre- 
dicted to occur in disordered systems. Also experiments 
in this spectral range can easily be done in the quan- 
tum limit, Tilo > ksT, with standard cryogenic tech- 
niques, i.e a pumped liquid A He optical cryostat (recall 
that 1 K ~ 20 GHz). 

The quasi-optical transmission measurement technique 
takes advantage of broadband sources of coherent electro- 
magnetic radiation known as Backward Wave Oscillators 
(BWO). The frequency of the output radiation from these 
sources can be swept by adjusting the applied potential. 
The resulting output power spectrum is extremely repro- 
ducible Recording the power received at the detector for 
1) the sample and then 2) a blank or free space, and 
taking their ratio provides the transmission data as a 
function of frequency. The striking feature of the trans- 
mission spectra is the periodically spaced peaks known as 
Fabry-Perot resonances which occur whenever the thick- 
ness of the sample is equal to an integer number of half 
wavelengths inside of the sample. The transmission can 
be recorded and the resonant peak's center frequencies 
and heights can be determined. From these measured 
parameters one can directly determine the real and imag- 
inary components of the complex conductivity. 

Below a certain temperature, it was found that 
the Fabry-Perot resonance transmission traces did not 
change. This necessarily means that below some temper- 
ature the real and imaginary components of the complex 
conductivity, become temperature independent. At low 
temperatures thermal contractions, which can affect both 



the peak heights and center frequencies become negligi- 
ble as well. Figure 5 shows this trend by plotting the 
conductivity of the 69% sample versus temperature for 
various constant frequencies with the DC conductivity. It 
is clear that at high temperatures the AC conductivity 
curves i ipproach and follow the temperature dependent 
DC conductivity, and the greater the frequency of a par- 
ticular curve, the greater the temperature at which this 
occurs. It is also clear that the iso-frequency curves do 
indeed flatten out at low temperatures, indicating that 
the T = K photon-assisted hopping behavior is being 
observed. All of the other samples showed a similar flat- 
tening out of the AC conductivity at low temperatures. 
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FIG. 5. Both the DC conductivity and the real component 
of the complex conductivity at various constant frequencies 
versus temperature on logarithmic scales for the 69% sample. 

Using DC data alone, evidence of the effects of 
electron-electron interactions, namely a crossover from 
Mott to ES type VRH is questionable at best. The range 
of temperature across which the various conduction pro- 
cesses might be occurring are minimal and it is difficult 
to have confidence that a fit to Mott-type 7 1-1 / 4 behavior 
or ES-type T~ x l 2 behavior is correct. The distinguishing 
features of an interacting versus a non-interacting disor- 
dered system are much easier to tell apart when inves- 
tigating the frequency dependent conductivity though. 
Specifically, linear and quadratic power law are clearly 
discernable and evidence for such a crossover in the fre- 
quency dependent conductivity spectrum is clearly seen 
as shown in Figures 6 and 7. 



IV. ANALYSIS AND RESULTS 

A. DC Conductivity 

Much of the story that the DC conductivity can tell 
is summarized in Figure 8 which shows the DC conduc- 
tivity of the 50% Si:P sample. The data in the upper 
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FIG. 6. The measured 2.8 K value of the real part of the 
complex frequency dependent conductivity for the Si:P 62% 
and 56% samples plotted versus frequency on logarithmic axes 
displaying a crossover in the type of conduction mechanism. 
The dashed line is a fit to the lower portion of the data and 
follows a nearly linear power law. The solid line through the 
upper portion of the data is a fit of that data and follows a 
nearly quadratic power law. 



FIG. 7. The measured 2.8 K value of the real part of the 
complex frequency dependent conductivity for the Si:P 45% 
and 39% samples plotted versus frequency on logarithmic axes 
displaying a crossover in the type of conduction mechanism. 
The dashed line is a fit to the lower portion of the data and 
follows a nearly linear power law. The solid line through the 
upper portion of the data is a fit of that data and follows a 
nearly quadratic power law. 



panel is plotted on a T -1 / 4 axis in order to assist vi- 
sual recognition of Mott type VRH behavior. The lower 
panel shows the sample data plotted on a T~ x / 2 axis in 
an attempt to elucidate any region of data that might 
be following ES type VRH behavior. When plotted on 
any inverse temperature scale, the left hand side repre- 
sents the high temperature region, and it is this higher 
temperature region, up to approximately 25 K, that has 
been fit to a linear function in the upper panel. How- 
ever, closer inspection of Figure 8 reveals that in the 
upper panel, although the higher temperature data has 
been fitted to T -1 / 4 behavior, the lower portion seems 
to be linear as well indicating a possible T -1 / 4 depen- 
dence as well on this semi- log plot. In the lower panel 
of Figure 8, the lower temperature data is fitted to a lin- 
ear function, possibly indicating that in this region, i.e. 
below T x « 12 K, the data might well be described as 
following ES type behavior. As in the upper panel, the 
other region of data, in this case the higher temperature 



data, shows a somewhat linear dependence at least over 
a finite range of temperature. 

That the upper and lower portions of the data seem 
linear on either scale summarizes the vicissitudes of try- 
ing to determine Mott or ES type behavior from the DC 
data alone. From Figure 8 one sees that the Si:P 50% 
sample has a crossover temperature of approximately 12 
K, whereupon the DC conductivity begins to follow a 
different dependence. It is unclear if this temperature 
sets the energy scale at which a crossover from Mott to 
ES type VRH occurs, but nevertheless it is clearly an 
important energy scale in this model disordered system. 

Further analysis shows that fitting the high tempera- 
ture and low temperature portions of Ode respectively to 
Mott and ES-type functional forms, across our experi- 
mentally accessible temperature range, is incorrect. The 
characteristic Mott temperatures, T , determined for our 
Si:P samples is of the order 10 7 K, a ridiculously large 
value. The ES characteristic temperatures, Tq, deter- 
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FIG. 8. The DC conductivity for the 50% Si:P sample plot- 
ted versus T" 1 / 4 in the upper panel and plotted versus T" 1 / 2 
in the lower panel. 

mined by fitting the lower temperature data are more 
reasonable but still extremely large, of the order 10 3 K. 
Compared to the characteristic Mott and ES tempera- 
tures measured by Dai et al. 26 and Hornung et al. 27 , our 
values seem extremely high and unphysical, even though 
their materials were of a different composition, or fell in 
a different range of dopant concentrations. 

Further evidence that the Mott and ES functional 
forms are inappropriate can be gleaned from analysis of 
the localization lengths from the characteristic temper- 
atures. Using the relationship for the ES characteris- 
tic temperature, i.e. Eq. 9 and the measured values 
of Tq, the localization lengths for the various Si:P sam- 
ples within our doping range come out to be less than 
2 nm, not an unreasonable length scale as the inter- 
atomic spacing in the diamond structure lattice of Si 
is 2.35 A, but an order of magnitude smaller than the 
20 nm found from analyzing the AC conductivity data, 
to be presented in the following section. With an ap- 
propriate estimate of the noninteracting density of states 
Nq(Ef), an estimate of the localization length can also 
be determined from Tq. For example the 62% sample, 
with a T = 39 x 10 6 K, estimating N (E F ) to be the 
ratio of the dopant atomic density x = 2.173 x 10 2 4 m" 3 



to an appropriate bandwidth, 0.1 eV, gives a localiza- 
tion length of approximately 6 A, a factor of 4 smaller 
than the one determined from Tq. That the localization 
lengths do not agree gives further evidence that 1) either 
our assumptions of what regions of the DC conductivity 
correspond to Mott or ES type VRH are incorrect or that 
2) deep in the insulating regime behavior different than 
that predicted for canonical Coulomb and Fermi glasses 
exists. 

The former case may well be true, (although evidence 
from the AC conductivity measurements seem to indicate 
otherwise), and these issues will be discussed presently, 
but first the latter issue shall be addressed. In fact 
such a prediction that deep in the insulating regime 
an alternate method of charge transfer should exist is 
exactly what Hornung et al. 27 claim to have observed 
in their DC conductivity measurements of Si:P. Among 
other findings, these authors claim that below a major- 
ity dopant concentration of 2.78 x 10 18 cm" 3 , (recall that 
x c = 3.52 x 10 18 cm" 3 ) the conductivity was better de- 
scribed by an activation process with an activation en- 
ergy E 2 which was interpreted as a "hard" Hubbard gap. 

Inspection ofthc 1.80xl0 18 cm" 3 and2.40xl0 18 cm" 3 
dopant concentration data presented in literature 27 , 
which were analyzed as an activation process, show that 
the data spanned a very small temperature range. Our 
69% sample corresponds to a dopant concentration very 
similar to the 2.40 x 10 18 cm" 3 dopant concentration 
sample in Ref. 29. Our DC conductivity results extended 
to 4.2 K. As indicated in Figure 3 a kink in the curve is 
clearly evident, which occurs at or somewhat below the 
lowest temperature indicated in the work from Hornung 
et al. 27 . The change in functional form that we observe 
is inconsistent with the presence of a hard gap which 
would result in a single activated process as the conduc- 
tion mechanism. As we will show below, further evidence 
for the absence of a hard gap comes from the fact that we 
observe appreciable AC conductivity at frequency scales 
well below that of the purported gap energy. 

B. AC Transport and Polarizability: The Crossover 

Some of the most stunning and germane results from 
these investigations of electron-electron interactions in 
disordered systems are best summarized visually in Fig- 
ure 9 in which the measured real part of the complex 
conductivity for the 69% and 50% samples are displayed 
as a function of frequency. The solid lines arc linear and 
quadratic fits to the lower frequency and higher frequency 
data respectively, qualitatively displaying the expected 
behavior for a crossover from a Coulomb glass to a Fermi 
glass. The overlayed dotted lines are fits to the data 
using Eq. (12) and clearly this formula provides only 
a rough guide. The fit was made by forcing the linear 
portion to pass through the origin as well as the low fre- 
quency data and leaving the prefactor of the quadratic 
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portion as a free parameter. The crossover between lin- 
ear and quadratic behavior is much more abrupt than 
the ES function predicts and is observed over our entire 
doping range. Such behavior was previously observed in 
an an analogous disordered semiconductor system, Si:B, 
for two very closely spaced doping levels near the MIT 3 . 
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FIG. 9. The measured real part of the complex conduc- 
tivity for the 69% and 50% Si:P samples as a function of 
frequency. The solid lines are linear and quadratic fits to the 
lower frequency and higher frequency data respectively quali- 
tatively displaying the expected behavior for a crossover from 
a Coulomb glass to a Fermi glass. The overlayed dotted lines 
are fits to the data using Eq. (12). 

The sharpness of the crossover is interesting in its own 
right and deserves further investigation as we do not here 
propose an explanation. This crossover in the real part of 
the complex conductivity was observed across our entire 
dopant concentration range except for samples very close 
the MIT, where the crossover presumably occurs below 
our experimental frequency range. The experimentally 
determined crossover frequency, lo Ci plotted versus con- 
centration in a functional form ^1 — -^j is summarized 
in Figure 10. 

The functional form in Figure 10, namely a power law 
in 1— — , is chosen because both the localization length, £, 
and the dielectric susceptibility of a sample near the MIT, 
namely 4nx = e\— est, are either predicted to or have 
been experimentally shown to follow such a form. Fitting 
the crossover frequency to such a functional form we find 
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FIG. 10. The experimentally determined crossover fre- 
quencies for the various Si:P samples. The data is plotted 
against concentration in a functional form (l — -^-)^ ■ We 
find P = 1.65. 



that lo c is proportional to \1 — where (3 = 1.65. By 

comparing the concentration dependencies of the relevant 
energy scales to this experimentally determined power 
we will show that the Coulomb interaction energy deter- 
mines the crossover. 

Recall that two energy scales have been set forth as 
being the energy scale at which the crossover from Fermi 
glass to Coulomb glass- like behavior is to occur. The 
original theoretical framework set forth by Efros and 
Shklovskii 2 has the Coulomb interaction energy U as the 
important energy scale. The second relevant energy scale 
in these systems is the Coulomb gap width, A. The evi- 
dence that the Coulomb gap is the pertinent energy scale 
of the system stems from an experiment performed by 
Lee et al. 3 . This report claims that it is not a gap in the 
single particle density of states, but rather a "dressed" 
Coulomb gap in the multi-particle density of states that 
determines the critical energy scale. Their claim, for 
which there does exist a theoretical framework 28 , is that 
hopping phenomena occur due to multi-particle processes 
or polaronic effects, but that a Coulomb gap still forms 
in the density of states of these charge carrying entities. 
Nevertheless it is not unreasonable to believe that the 
dressed Coulomb gap should take a similar form to that 
of the ES single particle gap, Equation (7), as both are 
presumably caused by the same long range Coulomb in- 
teraction. 

The dielectric constant enters into many of the ger- 
mane formulae describing the electrodynamics of these 
disordered insulating systems. The effects of screening 
and polarization become more and more relevant as the 
MIT is approached. Therefore it is of fundamental im- 
portance to develop a theory that self-consistently takes 
into account the effects of the imaginary component of 
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FIG. 11. The T = 2.8 K full measured dielectric constant 
ei for the various Si:P samples are plotted versus 1 — The 
inset is that portion of the dielectric constant due to the disor- 
dered interacting system itself, A-k\ obtained by subtracting 
off the background dielectric constant esi (of the host ma- 
terial silicon) from the full dielectric constant e\. A clear 
divergence as a function of concentration is seen. These mea- 
surements are consistent with previous reports, and facilitate 
the analysis of the crossover energy scale. 



the frequency dependent conductivity, i.e. the dielectric 
constant. As mentioned before, the immense advantage 
of our experimental techniques is the fact that we were 
able to separately measure both components of the com- 
plex conductivity. The full measured dielectric constant, 
as determined at T — 2.8 K, which is independent of fre- 
quency across our accessible frequency range (50 - 1000 
GHz), is displayed in Figure 11. The inset shows that 
portion of the dielectric constant due to the disordered in- 
teracting system itself, which we refer to as the dielectric 
susceptibility of the disordered system, 47rx, obtained by 
subtracting off the background dielectric constant esi (of 
the host material silicon) from the full dielectric constant 

Ei- 
Measurements of the concentration dependent diver- 
gence of the dielectric constant are not new but having 
independently determined the values for our samples at 
these frequency ranges allows for more accurate analysis 
of the crossover energy scales. As a measure of confi- 
dence, when compared to previous experimental data on 
the dielectric susceptibility portion, namely 47rx, of the 
full dielectric constant for Si:P 29 , which was found to fit a 

form (2^ — l)^, with £ = 1.15, we obtain from our data, 
when plotted and fitted in a similar fashion, a power 
also very close to unity, namely 0.96. However, mea- 
surements of the full dielectric constant, E\, reported in 
the literature on Si:B 30 , a disordered insulator analogous 

to Si:P, were fitted to a functional form of ^1 — -^j , 
and the power was found to be (' — —0.71. We ob- 



tain from our data, when plotted and fitted in a simi- 
lar fashion, a power for S\ oc ^1 — jf^j of £' = —0.38, 

which is almost a factor of 2 different. This is perhaps ac- 
counted for by the fact that the materials have different 
dopant material. One must be cautious when analyz- 
ing the full dielectric constant in this latter manner as 
is found sometimes in the literature because there is a 
constant additive term due to the background dielectric 
constant of the host material Silicon esi = 11-7 contained 
in the full dielectric constant E\ = Esi + 47rx, and only 
at dopant concentrations very close to critical does one 
obtain E\ « 47r% ;§> Esi- Our subsequent analysis of the 
crossover energy between the non-interacting Fermi glass 
behavior and the electron-electron interaction dependent 
behavior of the Coulomb glass state is carried out us- 
ing the latter normalized concentration functional form, 

— . It is therefore important to note here that as 

displayed in the inset of Figure 11, the dielectric suscepti- 
bility analyzed in this manner is found to obey the power 

law dependence, 47rx oc ^1 — jf^ , which presumably 

becomes the dominant power law dependence arbitrar- 
ily close to the MIT as the host material's esi becomes 
negligibly small. In fact such a trend can be seen in E\ 
in Figure 11, where at higher concentrations, namely the 
left-hand side of the figure, the slope tends to be steeper, 
as compared to the lower dopant concentrations on the 
right-hand side of the plot. 

Utilizing the measured values of the full dielectric con- 
stant for each sample one can now analyze the crossover 
energy. We mention briefly that in our analysis we pa- 
rameterize lo c in terms of concentration. We acknowledge 
that the crossover energy should not scale as a power law 
over the whole doping range, but in our limited range 
modelling shows this to be an acceptable parameteriza- 
tion. Now let us first examine the Coulomb gap width 
A. The experimental data shows that the crossover fre- 
quency lo c is proportional to ^1 — -^j where from our 
analysis (3 = 1.65, and the full dielectric constant E\ oc 
(1 — ^-1 . Using these values and Eq. 7 one finds 

that the DOS goes as No oc ^1 — jf^j . On the other 

hand the DOS can simply be estimated as No ~ x/W, i.e. 
the dopant concentration divided by the bandwidth W, 
which could only show the correct behavior stated above 
if the bandwidth were increasing faster than the dopant 
concentration itself. It is quite unlikely that the DOS 
could have a concentration dependence that would make 
this explanation consistent with the measurements, but 
one cannot entirely rule out that the resulting concentra- 
tion dependence is inconsistent with the Coulomb gap 
width as the DOS's concentration dependence is difficult 
to ascertain and the arguments provided above simply 
provide a reasonable estimate. 

There is strong support however that the crossover en- 
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ergy is determined by the Coulomb interaction energy 
though. By setting the measured crossover energy scale 
equal to the Coulomb interaction energy we are able to 
determine both the magnitude of the localization length 
and the exponent with which it diverges as a function of 
concentration. Using an appropriate pre-factor for the 
overlap integral 20 , Io = 10 13 s -1 , in the most probable 
hop distance term, r u = £[£n(2Io/?uv], we find a localiza- 
tion length dependence as shown in Figure 12. 
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FIG. 12. Calculated localization length plotted versus 
1 — — as determined from the measured crossover fre- 
quency energies. Our analysis entailed using an appropriate 
pre-factor for the overlap integral 20 , Iq = 10 13 s _1 in the most 
probable hop distance term, — (\ln{2Io/Uu>\. 

The localization length exponent as shown in the figure 
is close to unity, the value originally predicted by McMil- 
lan in his scaling theory of the MIT 10 , and the magnitude 
of the localization length is reasonable. The choice of the 
prcfactor for the overlap integral may vary, but due to the 
fact that it is located within a natural logarithm term, 
even an order of magnitude variation greater or smaller 
than 10 13 will result in a change of approximately a factor 
of two in the localization length. The consequent anal- 
ysis for the localization length exponent, namely after 
changing I by an order of magnitude both greater and 
smaller, still results in an exponent close to unity for both 
cases. It must be mentioned that some experimental re- 
sults have indicated a localization length exponent closer 
to 2) but these were found on the metallic side of the 
MIT using the zero temperature conductivity 29 . These 
results strongly point towards the Coulomb interaction 
energy as being the energy scale at which the observed 
frequency dependent crossover from ES to Mott-like hop- 
ping conduction occurs. 

Figure 13 offers a summary of various experimental 
data on different types of disordered insulating systems, 
including our own, of measured DC and AC crossover 
energies versus the normalized concentration parameter, 
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FIG. 13. Summary of various experimentally determined 
crossover energy scales determining the boundary between the 
interacting and non-interacting states of an electron glass sys- 
tem. The data represented by the open diamonds is adapted 
from tunnelling measurements 30 and the data represented by 
the full diamonds is adapted from AC conductivity network 
analyzer measurements performed in a dilution refrigerator 3 . 
Both measurements were done on Si:B, a disordered insulator 
analogous to Si:P. Our measured crossover frequencies 4 are 
represented by the full circles and our experimentally deter- 
mined crossover temperatures are converted to frequency and 
are represented by the crossed squares. 



^1 — -^j . The crossover energy for the various exper- 
imental techniques are converted to frequency and dis- 
played in the figure as GHz. The data in the upper left 
of the figure displayed as open diamonds is the Coulomb 
gap width in the single particle DOS as determined by 
tunnelling measurements in Si:B 30 , a disordered insulator 
analogous to Si:P. Although the Coulomb gap, speculated 
by Lee ct al. to set the scale of the linear to quadratic 
crossover, was a smaller re- normalized or "dressed" one, 
on this plot the single particle gap sets an overall scale. 
One can see that the crossover energy is well inside the 
single particle Coulomb gap for all energies. The data 
in the lower left portion of the figure depicted by the 
full diamonds are crossover frequencies as determined by 
Lee et al. from AC conductivity measurements using 
a network analyzer and a dilution refrigerator in order 
to achieve quantum limit-like behavior on the same set 
of Si:B samples as the above tunnelling data 3 . In their 
measurements, a sharp crossover analogous to our own 
observations was seen. The solid circles in Figure 13 
represent the crossover frequencies as determined from 
our frequency dependent conductivity measurements on 
Si:P. Our data extends the observed range of the sharp 
crossover deep into the insulating regime as the Si:B mea- 
surements were closer to the critical concentration. The 
Si:P and Si:B data do not fall on the same line. This may 
be due to slight differences in these various disordered 
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systems' dielectric constants or localization lengths. Fi- 
nally the open crossed squares represent the crossover en- 
ergies determined by us from our DC conductivity mea- 
surements. As can be seen the DC crossover data show 
a seemingly one-to-one correspondence to the observed 
crossover energies with the AC data. This in itself is a 
curious result, as the subsequent analysis of the DC data 
indicated that the observed crossover was not consistent 
with a crossover from Fermi glass Mott-like behavior to 
Coulomb glass ES-like behavior. 

As a final note on the crossover data, recall that in the 
literature (see Ref. 29) one finds experimental DC con- 
ductivity data analyzed to indicate that in Si:P below 
a certain dopant concentration, approximately x/x c = 
80%, a "hard" Hubbard gap opens up. Having observed 
appreciable AC conductivity at frequencies below the 
purported hard gap energy, seems to rule out the pos- 
sibility of a "hard" Hubbard gap, as the optical conduc- 
tivity would be expected to drop precipitously to zero if 
such a gap were present. In fact it seems that the picture 
set forth by Efros and Shklovskii describing a crossover 
in the real part of the complex conductivity from a linear 
to quadratic dependence on frequency due to the effects 
of electron-electron interactions, seems to hold from very 
close to the quantum critical regime to deep within the 
disordered insulating electron glass system at least qual- 
itatively if not quantitatively. From the expected form 
of U(x), i.e. Eq. 11, and the reasonable behavior of the 
extracted £, it seems likely that the crossover is governed 
by the Coulomb interaction strength of a resonant pair as 
predicted by Efros and Shklovskii, but the sharper than 
expected behavior still needs to be explained. 
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FIG. 14. The low energy, linear in frequency portion of 
the real part of the complex conductivity measured for three 
Si:P samples plotted versus frequency. This data focuses on 
the Coulomb glass portion of electrodynamic response. The 
dashed lines are forced linear fits indicating self-consistency 
with the theory for a Coulomb glass, namely a linear fre- 
quency dependence and an increase in the magnitude of ai 
with increasing dopant concentration. 



C. AC Transport and Polarizability: The Coulomb 
Glass 

Having established the fact that indeed a crossover in 
the real part of the frequency dependent complex con- 
ductivity occurs across a broad range of dopant concen- 
trations in an archetypical model system of a disordered 
insulator, we wish to focus on the low energy limit of 
this electron glass system, i.e. the Coulomb glass. We 
will attempt to highlight successes and shortcomings of 
the theoretical predictions by drawing on the results of 
both the real and imaginary components of the com- 
plex conductivity of Si:P and the previously reported 
results of another experimentally determined Coulomb 
glass, namely NbSi 5 . Comparing and contrasting the re- 
sults of two Coulomb glasses, an amorphous highly doped 
metal-semiconductor alloy, NbSi, and a comparatively 
sparsely doped crystalline semiconductor, Si:P, evidence 
for a general phase diagram of the Coulomb glass and 
the quantum critical region is revealed relative to our 
experimentally accessible window. 

Figures 14 and 15 show the low energy magnitudes of 
the real and imaginary components of the complex fre- 
quency dependent conductivity plotted against frequency 



for a number of Si:P samples. These figures focus on 
the low frequency Coulomb glass data from the measure- 
ments of our experimentally accessible region of the op- 
tical conductivity. The dashed lines are linear fits to the 
data, highlighting two telltale characteristics of Coulomb 
glass-like behavior, namely a linear dependence of both 
components of the complex conductivity on frequency, 
and an increase in the magnitudes of both a\ and o-i with 
increasing dopant concentration. The error bars shown 
in Figure 14 are representative of all the data point's er- 
ror bars; they show a general trend of increasing at lower 
frequencies both due to the logarithmic plot and the in- 
creasing difficulty in measuring such low conductivities. 
By contrast the error bars are much smaller and there is 
little scatter of the data points from the linear fits in Fig- 
ure 15 because due to the inherent nature of an insulating 
system, the relative magnitude of the imaginary compo- 
nent is greater than the real component of the complex 
conductivity, and therefore falls in a range that can be 
measured more accurately and precisely. 

Closer inspection of Figure 15 will reveal that the 
iT2 data is linear across the entire measured spectral 
range, whereas the previous section clearly established a 
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16. A model function of o\(lo), the dotted curve, 
and the Kramers- Kronig calculated magnitude of a<i{uj), the 
dashed curve, on axes of arbitrary conductivity and frequency. 
The model function for ai(u>) uses the Heaviside function to 
change from a linear dependence on frequency to a quadratic 
dependence at 10 on the arbitrary frequency scale. Above 100 
on the arbitrary frequency scale a Lorentzian function is used 
to simulate a peak in the optical conductivity spectrum. The 
shaded gray line indicates a linear dependence. 



FIG. 15. The low energy, linear in frequency magnitude 
of the imaginary part of the complex conductivity measured 
for four Si:P samples plotted versus frequency. This data 
focuses on the Coulomb glass portion of electrodynamic re- 
sponse. The dashed lines are forced linear fits indicating 
self-consistency with the theory for a Coulomb glass, namely 
a linear frequency dependence and an increase in the magni- 
tude of ct 2 with increasing dopant concentration. 

crossover in g\ as a function of frequency. One might 
question whether this is physically acceptable on the 
grounds of Kramers-Kronig (KK) compatibility. This is 
established in Figure 16. This figure shows model data 
for (Ti, plotted as the dotted line, on an arbitrary con- 
ductivity axis that via proper utilization of the Heavi- 
side function mimics the Fermi glass to Coulomb glass- 
like crossover, namely a change in the frequency depen- 
dence of the real part of the complex conductivity from 
approximately linear to approximately quadratic, at the 
value of 10 on the arbitrary frequency scale. Above 100 
on the arbitrary frequency scale a Lorentzian function 
is used to simulate a peak in the optical conductivity 
spectrum. The KK evaluated magnitude of cr 2 , result- 
ing from the model data of o\ , is shown with the dashed 
curve. The imaginary component clearly displays a lin- 
ear dependence on frequency, as evidenced by the shaded 
gray line, which is an aid to the eye for linear dependence, 
even though there exists an abrupt kink in the real part 
of the complex conductivity. 

Having established that below the crossover frequency, 
both components of the complex conductivity follow an 
approximately linear frequency dependence allows us to 



analyze the full complex conductivity, a, using the fol- 
lowing Kramers-Kronig compatible form, 

a = u\ + i(T2 = A(iu) a 
= Aio a cos (™) + iAu a sin (™) . (14) 

Such a dependence requires that the real and imaginary 
components follow the same power a over a broad range 
of frequency and that this simple KK-compatible form 
only holds for a < 1, thus this is perfectly valid when an- 
alyzing the Coulomb glass regime. Given the functional 
form of the complex conductivity in Equation (14), one 
can very accurately determine the power of a by taking 
the ratio of the magnitude of 02 and o\ to obtain what 
is commonly known as the phase angle of the complex 
conductivity 

a= l tan -ifN). (15) 

7T \ (71 / 

The ratio of \01\lo\ and the subsequent power of a was 
determined for each of the Si:P samples and the results 
are summarized in Figure 17. A similar analysis was 
performed on the insulating electron glass NbSi 5 , and the 
results from this Coulomb glass-like material are shown 
in the figure for comparison purposes. Panel (a) in Figure 
17 shows the ratio of the magnitude of the imaginary to 
the real component of the complex conductivity for both 
Si:P and NbSi. 

It must be noted that the magnitude of the imaginary 
component plotted in the figure, for both Si:P and NbSi, 
is that portion of the polarizability attributed to the in- 
teracting disordered electron system itself. Recall that 
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FIG. 17. Panel (a) shows the ratio of the magnitude of <72 to 
<T\ for the two Coulomb glasses measured in this work, namely 
Si:P and amorphous NbSi. Panel (b) shows the calculated 
powers of a as determined from Equation (15). The dashed 
line through the NbSi data is a guide to the eye showing an 
approach to the value |, as expected upon approaching the 
quantum critical regime. Panel (c) shows the divergence of 
the magnitude of a\ as a function of concentration. 

the experimentally determined quantity is the full dielec- 
tric constant, E\ = £5^ + 4^%. The imaginary component 
of the complex conductivity plotted in the figure has been 
calculated from the dielectric susceptibility portion At:\, 
resulting from the electron-electron interactions via the 
relationship cr 2 °c (47rx) uj. The ratio of \<t 2 \/<ti us- 
ing an imaginary component of the complex conductiv- 
ity calculated from the full measured dielectric constant 
would necessarily diverge as the dopant concentration de- 
creased, i.e. deep in the insulating regime, because |cr 2 1 
would remain constant due to the static background di- 
electric constant from silicon itself, while o\ would con- 
tinuously decrease. 

Returning to panel (a) in the figure, we first note that 
this ratio for Si:P remains large and approximately con- 
stant across our range of dopant concentrations. From 
theory, one expects |<7i| to be approximately equal to 
| (T2 1 to within a factor of 2-5 (with a reasonable estimate 
for J ) as predicted by Efros 31 . Applied to Si:P, Eq. (12) 



in the T%uj < U(r w ) limit correctly predicts a linear cor- 
respondence between a\ and | cr 2 1 , but the theory incor- 
rectly predicts the measured proportionality by at least 
a factor of thirty. The proportionality is closer for NbSi 
but has a dependence on the doping concentration. This 
may be related to entering the quantum critical (QC) 
regime as discussed below. The magnitude of the ratio 
becomes even larger and the discrepancy greater if the 
full dielectric constant e\ is considered instead. 

Panel (b) in Figure 17 shows the power a as determined 
by Eq. (15) for both the amorphous NbSi and Si:P sam- 
ples. The values for Si:P are approximately equal to, but 
slightly less than one, consistent with Figures 14 and 15. 
This indicates that the prefactor of the real and imag- 
inary components of the complex conductivity for our 
samples of Si:P have the same concentration dependence 
across our entire range of dopant concentrations. 

Panel (c) in Figure 17 shows the magnitude of the 
real part of the complex conductivity for the Si:P sam- 
ples as the MIT is approached from the insulating side. 
This demonstrates that the prefactor A as set forth in 
Equation (14) can be written as a function of the nor- 
malized concentration, i.e. A oc (1 — — ) for Si:P. The 
prefactor of the imaginary component can be analyzed 
in a similar fashion and shows the same functional di- 
vergence as the real component, i.e. ~ (1 — :f-)~ 2 ' 3 - 
This indicates that the real and imaginary components 
have the same concentration dependence consistent with 
the simple Kramers-Kronig form of the complex conduc- 
tivity given in Equation (14). The fact that across a 
broad dopant range on the insulating side of the MIT we 
observe a similar divergence in both components of the 
complex conductivity for Si:P is in accordance with the 
concentration dependent theoretical predictions for pho- 
ton assisted hopping conductivity 2 , for excitations both 
within and larger than the Coulomb gap energy scale. 

Still focussing on the Coulomb glass regime, the exper- 
imentally determined concentration dependence of the 
prefactor A, namely (1 — j-) -2 ' 3 , as shown in panel (c) 
of Figure 17, can be used to determine which theoreti- 
cally predicted functional form for the real part of the 
complex conductivity is more germane for our measured 
system. Recall that for photon energies smaller than 
the Coulomb gap width, i.e. excitations outside of the 
gap, the concentration dependent terms are summarized 

by the following: ai(w) = Aui and that 

for the alternate case, hopping conduction inside of the 
Coulomb gap one finds: = Auj oc e\uj. Having ex- 

perimentally determined a parameterization for the con- 
centration dependencies for both the full dielectric con- 
stant, £1 oc (1 — :f-)~ ' 4 , and the localization length, 

£ oc (1 3L)-o.83 we g n j tnat tne f unc tional form for 

the prefactor is better described by photon-assisted hop- 
ping conduction occurring outside of the Coulomb gap. 
We find that £ 4 / £ i oc (1 — :f-)~ 2 ' 9 , which is not in exact 
agreement with our experimental findings for the prefac- 
tor (and in fact the density of states, No term would make 
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this a slightly steeper function), but this is in much better 
agreement than the predicted concentration dependence 
of the prefactor for photon assisted hopping conduction 
inside the Coulomb gap. 
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FIG. 18. A schematic of a possible phase diagram indicat- 
ing the relative positions of the crossover energy scale to QC 
behavior vs. concentration of Si:P and NbSi. The solid box 
indicates the "window" from which our experimental results 
are derived. 



Figure 18. The differently shaded regions indicate the rel- 
ative positions of the crossover energy scale to QC behav- 
ior for Si:P and NbSi. As mentioned above the crossover 
between the QC and the Coulomb glass regimes occurs 
when the localization length scale becomes comparable 
to the dephasing length scale, namely £ « i w , and from 
this relationship one finds that the crossover frequency is 
inversely proportional to the DOS 9 . This crossover then 
is expected to be narrower in Si:P as compared to NbSi 
due to its smaller DOS. The crossover condition for the 
frequency in terms of the normalized concentration is 

u Q c*tl(l- (16) 

where the subscript QC refers to the quantum critical 
regime. The functional dependence is expected to be 
universal but the prefactor can vary from system to sys- 
tem. This then defines a condition for the crossover from 
QC to insulating glass-like behavior that depends only on 
the dynamical scaling exponent, z, and the localization 
length exponent, v, that is consistent with experimental 
evidence. For the case of z « 2, and v « 1 as evidenced in 
NbSi 6 and Si:P 4 respectively, the boundary would then 
have an approximately quadratic dependence on the nor- 
malized concentration as shown in Figure 18. 

V. DISCUSSION 



D. Quantum Critical Regime 

When approaching the MIT from the insulating side, 
the frequency dependence is expected to cross over to the 
QC behavior 6,32 , i.e. a\ oc w 1 / 2 in NbSi, when the local- 
ization length £ becomes comparable to the characteristic 
frequency dependent length scale, the dephasing length, 
i^ 12 , as discussed in Section II A. The crossover is not a 
phase transition and is not expected to be sharp. There- 
fore looking at a fixed window of frequencies, a broad, 
smooth crossover from u — > u 1 / 2 would show an aver- 
aged power of the frequency dependence similar to that 
measured for NbSi shown in the middle panel of Figure 
17. The fact that we see ana w 1 across our whole doping 
range in Si:P, but an a that approaches | in NbSi indi- 
cates that the critical regime in Si:P is much narrower. 
Simple dimensional arguments 33 give a result similar to 
the non-interacting case 9 that the crossover should be 
inversely proportional to the dopant DOS. The vastly 
smaller (a factor of 10 3 ) dopant density in Si:P relative 
to NbSi is consistent with a narrower QC regime in Si:P 
as compared to NbSi. 

A schematic outlining our experimental window and 
the boundaries for the crossover from the insulating 
Coulomb glass regime to the QC regime in both Si:P 
and NbSi versus a normalized concentration is shown in 



In analogy to the three phases of matter we are accus- 
tomed to dealing with in elementary physics, the solid, 
liquid and gas phases, a collection or assembly of elec- 
trons can be classified in a similar manner depending 
on the electrodynamic response of the system, each pos- 
sessing an accompanying theoretical model. The room 
temperature response of a metal upon applying an elec- 
tric field is well described by the Drude model. Here the 
electron system can be thought of as an electron gas with 
a temperature dependent mean free path and scattering 
time. An electron liquid would naturally be a system of 
electrons obeying Landau's theory of Fermi liquids 34,35 . 
This successful theory has distinct hallmarks in the low 
energy limit. We shall investigate this low energy behav- 
ior more closely in the following section by comparing it 
to the low energy behavior of the materials we studied, 
generally classified as electron glasses. 

The poorly understood and infrequently studied elec- 
tron glasses might be considered a genre of the final phase 
mentioned above, namely an electron solid, where the 
electronic states are localized due to disorder. The amor- 
phous nature of the glassy state connotes disorder and 
hence this appellation was chosen to describe this type 
of assembly of electrons. Via our unique and versatile AC 
conductivity measurement techniques, and by investigat- 
ing a series of disordered insulating systems proximal to 
a MIT, we have sufficient experimental evidence for elec- 
tron glass systems to propose a classification scheme, i.e. 
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a 'taxonomy' for the electrodynamic response of electron 
glasses, along with a phase diagram containing appropri- 
ate crossover boundaries for the electron glass, Coulomb 
glass and quantum critical regimes. 

A. On the Nature of the Coulomb Glass 

One of the fundamental reasons for undertaking our 
investigations of these electron glasses is the surprising 
fact that even in a material as thoroughly studied as 
doped bulk silicon (due to the relevance of such research 
for the semiconductor industry), there exists no clear 
consensus as to the ground state and the nature of the 
low energy excitations in these disordered insulating sys- 
tems. Using the frequency dependent complex conduc- 
tivity we wished to probe the question as to whether or 
not the Coulomb interactions fundamentally change the 
ground state of the system. A philosophy that perse- 
vered for many years was that in analogy to the Fermi 
liquid theory, the theory describing the excitations at ar- 
bitrarily low energies in a system localized due to dis- 
order, would indeed not be qualitatively influenced by 
electron-electron interactions 36 . The theory created un- 
der this pretense was the one set forth by Mott 1 , which 
gives an approximately quadratic frequency dependence 
of the real component of the complex conductivity. The 
coining of the appellation Fermi glass to describe this 
non-interacting disordered electron system was a direct 
consequence of the belief that there existed an analogy 
between the low energy electrodynamics of a Fermi liquid 
and this Fermi glass. Later, theoretical work and subse- 
quent investigations of the DC conductivity lent credence 
to the import of long range interactions stemming from 
negligible screening of the low energy excitations, e.g. the 
work by Pollak 14,28 and by Efros and Shklovskii 2 . 

In 'many' systems in solid state physics, this claim cer- 
tainly is true, namely that the effects of correlations be- 
come simpler at low frequency or energy. The canoni- 
cal example, as stated, is the Fermi liquid. The heavy- 
electron compounds such as CeAl% and UPtz are mate- 
rials that have strongly correlated, yet delocalized charge 
carriers, and in the low energy limit, they follow the non- 
interacting Fermi liquid theory with renormalized param- 
eters such as the effective mass and the relaxation rate 37 . 
For these heavy- fermion systems there is a one to one cor- 
respondence between the low-lying excitations and those 
in a nearly free electron gas when the parameters of the 
latter, namely the mass and scattering rate, are properly 
renormalized to take into account the correlations. The 
name 'heavy-electron' stems from this effective enhanced 
mass (often several hundred times) which is a direct con- 
sequence of the electron-electron interactions, but never- 
theless the theory describing the electrodynamic response 
of these systems is a non-interacting one. 

On the other hand, our experimental data for both 
NbSi and Si:P, two variants of electron glass systems, 



show Coulomb glass-like behavior, i.e. a theory depen- 
dent on interactions, in the low energy frequency depen- 
dent electrodynamic response. Experiments done by an- 
other group on an electron glass system 3 , namely Si:B, 
also show a linear in frequency power law dependence 
of the T = 0, uj — > conductivity supporting and veri- 
fying the import of electron-electron interactions in the 
low-energy response of these disordered systems. While 
this work on Si:B reported on only two concentrations 
of dopant very close to the MIT, our further investiga- 
tion across a broad range of dopant concentrations in 
Si:P can confirm that this linear power law dependence 
continues deep into the insulating regime. Thus we have 
clear evidence that an interacting theory, namely that of 
a Coulomb glass, dependent on electron-electron corre- 
lations, is correct at the lowest experimentally accessible 
energy scales in these electron glass systems, in opposi- 
tion to the original thoughts on the matter. Thus the 
Coulomb glass is a fundamental example in solid state 
physics as the non-interacting functional form is not re- 
covered in the asymptotically low energy limit. 

B. The Taxonomy of an Electron Glass 

Now we focus on the electron glass phase as a whole 
with the intent of setting forth a universal framework for 
the AC response of such disordered insulating systems. 
In summary, we have observed Coulomb glass-like be- 
havior across our entire range of doping concentrations 
in Si:P and via a comparison with NbSi have arrived at 
a boundary condition consistent with our experimental 
data for the crossover from QC to Coulomb glass-like be- 
havior, namely Equation (16). We have also established 
that a crossover from Fermi glass-like to Coulomb glass- 
like behavior exists, as is manifest in the experimental 
data of the real part of the complex frequency depen- 
dent conductivity, not only close to the MIT 3 , but also 
deep into the insulating electron glass regime. Employing 
these phenomenological results, namely that the power 
laws seen in the frequency dependent complex conduc- 
tivity are qualitatively consistent with the three expected 
types of frequency dependent behavior theoretically pos- 
tulated for the electron glass regime, we propose that 
the real component of the frequency dependent complex 
conductivity of these disordered insulators can be param- 
eterized in general by 

en oc tu a ^^ (17) 

where the power a depends on the frequency and the 
dopant concentration. This classification scheme, or tax- 
onomy, based on a(u),x), captures the essential and ger- 
mane physics of the electrodynamic response of electron 
glasses. 

Figure 19 summarizes our proposed phase diagram and 
taxonomy on a plot of frequency versus the normalized 
concentration (1 — ^-). Recall that the criterion for 
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VI. CONCLUSION 
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FIG. 19. This schematic of a possible phase diagram for 
an electron glass shows the regime boundaries for the QC, 
the Coulomb glass and the Fermi glass regimes in a two di- 
mensional space of frequency versus normalized concentration 

the regime boundary condition for the crossover from 
Coulomb glass or Fermi glass-like behavior to the region 
where QC dynamics dominated was determined by Eq. 
(16). For a z = 2 theory and a localization length expo- 
nent v = 1, this crossover would go as uj qc oc (1 — -^-) 2 
with a prefactor dependent on the dopant DOS. fbr 
the crossover from Fermi glass to Coulomb glass-like be- 
havior, shown on the right hand side of the figure as 
the light gray and white regions respectively, we found 
that the crossover energy was more likely to be deter- 
mined by the Coulomb interaction energy as opposed to 
the Coulomb gap width using concentration dependent 
arguments. Either way our experimentally determined 
crossover had the following concentration dependence, 
uj c oc (1 — j-) 1 ' 65 - Here we must note that though the 
curve depicting this latter crossover, as shown in the fig- 
ure, might intersect the QC crossover curve (depending 
on the material's energy scales), the QC regime will al- 
ways take precedence. Recall that this crossover energy 
stems from a length scale argument that quantum critical 
behavior occurs when the electronic localization length, 
£, is greater than the dephasing length, £ u . Ergo since 
both the Fermi glass and Coulomb glass behavior only 
occur in the opposite case, e.g. £ < £ u , the length scale 
dependent crossover condition for QC behavior always 
dictates the low energy, close to critical behavior. 

The classification of the power of o\ oc u> a ^' x \ e.g. 
a = 1/2 for QC dynamics, a = 1 for Coulomb glass be- 
havior and a = 2 for Fermi glass behavior, provides more 
than just a taxonomy for these disordered insulating sys- 
tems. For the latter two regimes, with powers 1 and 2 
respectively, these numbers also represent the dimensions 
of phase space changing directly as Uui that these systems 
have access to for the creation of a real excitation. 



In summary, we have presented results of our frequency 
dependent complex conductivity investigations into cor- 
related electron phenomena in the electrodynamics of the 
electron glass Si:P. The insulating side of the MIT dis- 
played the expected power law dependencies for a Fermi 
glass and a Coulomb glass system with the latter be- 
ing the low energy limit in our experimentally accessi- 
ble window. In comparison with another electron glass, 
NbSi, with a distinctly different DOS, wherein we see a 
crossover from Coulomb glass-like behavior to quantum 
critical behavior, we have formulated a phase diagram 
based on phenomcnological observations of all the pre- 
dicted states for the insulating side of the MIT. The tax- 
onomy for the electron glass that we have presented, is 
an attempt to both capture the essence of the physics 
and allow for a general description of the low energy fre- 
quency dependent electrodynamics of electron glasses. 

A number of issues have also been left unresolved by 
this work. First of all, the apparent sharpness of the 
crossover in the real part of the frequency dependent con- 
ductivity cannot be explained by the current model of 
Efros and Shklovskii 2 . Second, there certainly also exist 
further misgivings with respect to the theories concerning 
the DC conductivity. For instance, the DC conductivity 
data does not seem to fit the Mott or ES predictions par- 
ticularly well for our doping levels well into the insulating 
phase as indicated by the unphysical characteristic tem- 
peratures derived therefrom. Also, we have experimen- 
tally observed a crossover energy in the DC conductivity 
that is quite comparable to the crossover energy seen 
in the frequency dependent conductivity, and there is no 
reason within the canonical theory that there should be a 
correspondence in the AC and DC energy scales. Thirdly 
the experimentally determined ratio of the imaginary to 
the real component of the complex conductivity shows 
a deviation from the predicted ratio. Should a subse- 
quent re-analysis of the theories pertaining to the elec- 
trodynamics of these electron glasses be undertaken, one 
should certainly include a self-consistent dynamic dielec- 
tric constant, and in a correct description, the proper 
ratios and relative sharpness of the crossovers must be 
manifest. 

We are also hopeful that this work will plant the seed 
for further future investigations of other disordered insu- 
lating systems. For instance the samples used in these in- 
vestigations were nominally uncompensated. Purposeful 
compensation will naturally redistribute many relevant 
energy scales. It would be of fundamental importance to 
investigate the veracity of the proposed taxonomy of the 
electron glass with respect to compensated systems. One 
must also note that our phase diagram stemmed from 
phenomenological observations of the three types of elec- 
tron glass behavior, namely Fermi glass, Coulomb glass 
and quantum critical, but we would be remiss if we ne- 
glected to emphasize that we did not physically observe a 
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crossover from Fermi glass to quantum critical behavior, 
just the other crossovers that are possible in the taxon- 
omy phase diagram. Needless to say, such an observation 
would be of great import in fully understanding electron 
glass behavior. Finally we have seen that in our experi- 
mental window, the original expectations set forth in the 
literature 36 have proven to be not realized, in that at our 
lowest accessible frequencies, an interacting theory is the 
germane one. This does not however rule out a subse- 
quent change in the conduction mechanism at even lower 
frequencies, such as non-linear effects. 



VII. ACKNOWLEDGEMENTS 

We would also like to thank S. Kivelson, B. Shklovskii 
and E. Abrahams for useful discussions and also N. F. 
Mott who introduced one of us (GG) to the subject. This 
research was supported by the National Science Founda- 
tion grant DMR-0 102405. 



25 G. Griiner, Millimeter and Subrnillimeter Wave Spec- 
troscopy of Solids (Springer Verlag, Berlin, 1998). 

26 P. Dai et al., Phys. Rev. Lett. 66, 1914 (1991). 

27 M. Hornung et al., Physica Status Solidi (b) 218, 75 (2000). 

28 M. Pollak and M. Ortuiio, in Electron-Electron Interactions 
in Disordered Systems, edited by A. L. Efros and M. Pollak 
(Elsevier, New York, 1985), pp. 287-408. 

29 H. F. Hess et al., Phys. Rev. B 25, 5578 (1982). 

30 M. Lee et al., Phys. Rev. B 60, 1582 (1999). 

31 A. L. Efros, Sov. Phys. JETP 62, 1057 (1985). An 
analagous three-dimensional form of the two-dimensional 
theory described in this paper can be derived. 

32 H.-L. Lee et al., Phys. Rev. Lett. 80, 4261 (1998). 

33 S. A. Kivelson, private communication. 

34 L. Landau, Sov. Phys. JETP 3, 920 (1957). 

35 L. Landau, Sov. Phys. JETP 8, 70 (1959). 

36 P. W. Anderson, Comments S. S. Phys. 2, 193 (1970). 

37 A.M. Awasthi et al., PRB 48, 10692 (1993). 



1 N. F. Mott and E. A. Davis, Electronic Processes in Non- 
Crystalline Materials, Second Edition, (Oxford University 
Press, Oxford, 1979). 

2 A. L. Efros and B. I. Shklovskii, in Electron- electron Inter- 
actions in Disordered Systems, edited by A. L. Efros and 
M. Pollak (Elsevier New York, 1985), p. 409-482. 

3 M. Lee and M. L. Stutzmann, Phys. Rev. Lett. 30, 056402 
(2001). 

4 E. Helgren et al., Phys. Rev. Lett. 89, 246601 (2002). 

5 E. Helgren et al., Phys. Rev. Lett. 87, 116602 (2001). 

6 H.-L. Lee et al., Science 287, 633 (2000). 

7 F. J. Wegner, Zeitschrift far Physik B 25, 327 (1976). 

8 E. Abrahams et al., Phys. Rev. Lett. 42, 673 (1979). 

9 B. Shapiro and E. Abrahams, Phys. Rev. B 24, 4889 
(1981). 

10 W. L. McMillan, Phys. Rev. B 24, 2739 (1981). 

11 D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 
(1994). 

12 S. L. Sondhi et al., Rev. Mod. Phys. 69, 315 (1997). 

13 N. F. Mott, J. of Non-Cryst. Sol. 1, 1 (1968). 

14 M. Pollak, Disc. Far. Soc. 50, 13 (1970). 

15 V. Ambegaokar et al., Phys. Rev. B 4, 2612 (1971). 

16 A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 

17 V. L. Nguyen et al., Sov. Phys., Semicond. 13, 1281 (1979). 

18 A. L. Efros et al., S. S. Comm. 32, 851 (1979). 

19 S. Tanaka and H. Y. Fan, Phys. Rev. 132, 1516 (1963). 

20 Personal correspondence with B. I. Shklovskii. The Io used 
is determined by the Bohr energy of phosphorous. 

21 W. Teizer et al., Phys. Rev. Lett. 85, 848 (2000). 

22 W. R. Thurber et al., J. Elcctrochem. Soc. 127,1807 
(1980). 

23 Y. Ootuka et al., S. S. Comm. 36, 827 (1980). 

24 A. Schwartz et al., Rev. Sci. Instrum. 66, 2943 (1995). 



18 



